DIFFERENTIATION EVENS OUT ZERO SPACINGS 



DAVID W. FARMER AND ROBERT C. RHOADES 

Abstract. If / is a polynomial with all of its roots on the real line, then the roots of the 
derivative /' are more evenly spaced than the roots of /. The same holds for a real entire 
function of order 1 with all its zeros on a line. In particular, we show that if / is entire of 
order 1 and has sufficient regularity in its zero spacing, then under repeated differentiation 
the function approaches, after normalization, the cosine function. We also study polynomials 
with all their zeros on a circle, and we find a close analogy between the two situations. This 
sheds light on the spacing between zeros of the Riemann zeta-function and its connection 
to random matrix polynomials. 



1. Introduction 

If / is a polynomial with all of its zeros on the real line, then all of the zeros of the 
derivative /' also lie on the real line. The same holds for entire functions of order 1 which 
are real on the real line. But while the zeros of / in an interval could in principle be any 
finite subset of the interval, the condition that / have only real zeros will impose restrictions 
on the zeros of /'. The main idea of this paper is that if the zeros of / lie on a line then the 
zeros of /' are more evenly spaced than the zeros of /. Repeated differentiation leads the 
zeros to become more and more evenly spaced, and under appropriate conditions the zeros 
will approach equal spacing. See Section 1231 for precise statements. 

There is a long history to the behavior of the zeros of entire functions under repeated 
differentiation. Polya ^B] made several interesting conjectures for the case of real entire 
functions, some of which have only been solved recently El El El El- Many of 
those questions concern the reality of zeros of a real entire function, or the motion of zeros 
toward (or away from) the real line when the function is differentiated. In this paper we 
take a complementary view and consider the motion along the real line when the function is 
differentiated. 



1.1. Zeros of the derivative. It is tempting to think of the zeros of /' as lying close to the 
midpoint of neighboring zeros of /. That is a useful model, and we analyze the sequence of 
midpoints and use it as a tool to study the zeros of the derivative. A more accurate picture is 
think of the zeros of /' as trying to move as far as possible from the zeros of /, and so move 
toward regions where there are fewer zeros of /. See Figure ^ for an example. If Zj < zjj^i 
are consecutive zeros of /, then RoUe's theorem asserts that there will be a zero z' of /' 
between zj and ^j+i- If there are a large number of other zeros near zjj^i, then z' will be 
closer to zj. 
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There are several intuitive reasons for this phenomenon. Suppose Zj are the zeros of /, 
listed in increasing order, and consider 

(1.1.1) jW-^T^- 

If Zj < p < Zj+i then f'{p) = if and only if ^(p) = 0. Substituting z = p into we 
see two kinds of terms: positive for k < j, and negative for k > j + 1. If there are many 
zeros near Zj^i then all of those negative terms have to be balanced by something. If there 
aren't many zeros near Zj, then the only option is moving p closer to Zj. One can think 
of M = + Zj+i)) as a measure of how far the zero of /' will be from the midpoint. 

Specifically, 



1.1.2) 




Zj + Zj+i) + ^{zj^i ~ ZjY if M is small 
or Zj+i) — jj if M is large. 



An illustrative example is the polynomial / with zeros at {0, 1, ... , A^}. The first zero of 
the derivative is very close to 1/logA^, and (/'//)(!) ~ — logA^. 

A less explicit version of the above explanation uses Gauss' electrostatic model. By 1.111 . 
the zeros of /' can be thought of as points of zero electric field for a collection of equal charges 
located at {z^}, where like charges repel according to an inverse linear law. If there is a region 
of high concentration of charges, then these will push the equilibrium points toward regions 
of lower charge concentration. 

Another approach concerns the relationship between the spacing of zeros and the size of 
the relative maxima and minima of /. Roughly speaking, / gets big when there are larger 
than average sized gaps between zeros. But since / is a polynomial with all zeros on a line, it 
can't change slope too fast. This forces the locations of the relative maxima to shift around 
so that the slope of / can change appropriately. See Figure ^ particularly near the first local 
minimum and the third local maximum. 




Figure 1. The circles are the midpoints of neighboring zeros of the polynomial /, 
and the dots are the zeros of /'. Think of the circles as a target, which the zeros 
of /' often miss. 
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This approach is interesting because the current view is that the spacing of zeros is key 
to understanding the large values of functions like the Riemann zeta-function. 

In addition to entire functions with zeros on a line, we also consider polynomials with all 
their zeros on a circle. We find a close analogy between the two cases. The motivation here 
is that much recent work has concerned the connection between L-functions (such as the 
Riemann (^-function, which conjecturally have all their zeros on a line) and the characteristic 
polynomials of random unitary matrices (which have all their zeros on the unit circle). The 
connection between random matrices and L-functions is discussed extensively in Pj. 

In the next section we start with an old result that motivates the remainder of the paper. 
We first consider polynomials, and then generalize to entire functions of order 1. Our main 
results are stated in Section 12.41 In Section |21 we consider a somewhat simpler averaging 
process which we later will compare to differentiation. In Section 0] we consider the case of 
zeros which are almost equally spaced, and determine the rate at which differentiation evens 
out the spacing. Note that in Sections El and El we treat in parallel the cases of zeros on a 
line and zeros on a circle. In Section El we complete the proof of our main results, and also 
discuss the case that zeros do not lie on a line. 

We thank Robin Chapman, Brian Conrey, Nathan Feldman, Chris Hughes, and Paul 
McGuire for helpful conversations, and the referee for several useful suggestions. 

2. Small gaps become larger 

If / is a function with all its zeros on a line, define the "smallest gap between zeros" 
function 

G[f] = inf \aj - a^l 

where aj and are distinct zeros of /, with zeros repeated according to their multiplicity 
(so G[f] = if / has a multiple zero). We will show that if / is a polynomial with all zeros 
on a line, then G[f'] > G[f]. That is, the smallest gap between zeros of /' is at least as large 
as the smallest gap between zeros of /. This is in keeping with the underlying principle that 
zeros of /' tend to "move away" from regions with higher density of zeros of /, and to "move 
toward" regions with lower density of zeros of /. 

2.1. Motivating result. The following theorem is attributed to M. Riesz, see PHj- The 
result was rediscovered several times 1^ 1^ . including some results which are quanti- 
tatively stronger than the statement below. Our proof is similar in spirit, but somewhat 
simpler, than previous versions. 

Theorem 2.1.1. Let f{x) = e°'^h{x), where h is a polynomial with only real zeros. If a and 
a are real then G[f' + af] > G[f], with strict inequality provided the zeros of f are simple. 

In particular, zeros of /' are 'further apart' than zeros of /. 

It is implicit, and we give Lagrange's proof following equation ()2.1.2p below, that all zeros 
of /' + af are real. 

Proof. We may assume that / has only simple zeros qq < ai < ■ ■ ■ < aj . We can write 

J 

(2.1.1) /(x) = e"^n(^-a,) 

j=0 
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and so 

(2.1.2) ^^^{x) = a + a + y ^— . 

We see that all the zeros of /' + af are real because if x has positive (negative) imaginary 
part, then each term in the sum has negative (positive) imaginary part. Also, between each 
pair of zeros of / there is exactly one zero of /' + af, and there is one additional zero if 
a + a 7^ 0. 

Let q < p he successive zeros of /' + af. Assume, for a contradiction, p — q < Oj+i — aj 
for all j. Note that this is equivalent to 

1 1 

(2.1.3) > 



p — flj+i q — aj^ 



because p — aj+i and q — aj are either both positive or both negative, because there is only 
one zero of / between p and q. 



We have 



= a + a + y 

^-^ Tl — , 



— ' p — ai 



and 



J 1 
a + a +y 



j=0 



q-a^ 



Subtracting the two equations gives 

J 



j=0 



p — aj q — aj 



1 1 ■''^ 

(2.1.4) = — 

p- ao q- aj ~^ \P - ctj+i q - a^ 

The first two terms above are strictly positive, and by ()2.1.3j) all the terms in the sum are 
nonnegative. So the right side is strictly positive, which is a contradiction. □ 

Repeatedly applying the result leads to various linear combinations of / and its derivatives 
which satisfy the corresponding inequality. For example, G[f — cf"] > G[f] provided c > 0. 

2.2. Entire functions of order 1. We wish to generalize the above result to functions 
having infinitely many zeros. The appropriate class to consider is entire functions of order 
at most 1. We recall the standard terminology and properties of these functions. For more 
details, see jTj. 

If / is an entire function, then / is of order p if 

/oon r log log M(r) 
(2.2.1) limsup = p, 

r^oo log r 

where M(r) = max\z\=r \f{z)\. That is, f{z) = 0(exp(|z|''+'')) for all 5 > and no S < 0. 
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If / is of order p, then / has type r if 

(2.2.2) lim sup r"'' log M(r) = r. 

r— >oo 

That is, f{z) = 0(exp((r + 5)|2;|^)) for all 5 > and no 5 < 0. We say that / has minimal 
(respectively, finite or maximal) type if r = (respectively, 0<r<ooorr = oo). 

For example, cos(z) = |(e*^ + e~*^), so cos(z) has order 1 and type 1. And by Stirling's 
formula, l/r{z) has order 1 and maximal type. 

Of particular interest is the connection between order, type, and the distribution of zeros. 
Let (zn) denote the zeros of /, with Zj < zj+i and zeros repeated according to their multi- 
plicity. Denote by n{r) the number of zeros of / in \z\ < r. If / has only real zeros, then 
let n+(r) (respectively, rz_(r)) denote the number of zeros in (0,r] (respectively, [— r, 0)). A 
characterization of functions of finite type is due to Lindelof: 

Proposition 2.2.1. (Lindelof) If p is an positive integer, then the entire function f of order 
p is of finite type if and only if n{r) = 0{r^) and the sum over zeros 

(2.2.3) S{r) = J2 ^n" 

\Zn\<r 

is bounded. 

That cos(2;) has finite type and l/r{z) has maximal type follows directly from Lindelof 's 
theorem. 

The Riemann ^-function is not entire because it has a simple pole at s = 1, but (s — l)^(s) 
is entire. From the Dirichlet series representation ({s) = X^^i ^"'^ follows that the (- 
function is bounded for 3?(s) > 1 + 6. By the functional equation 

C(s)=2%^-isin(is7r)r(l-s)C(l-s) 

and the fact that the sine and F-function have order 1, it follows that the ^-function has 
order 1. 

The Riemann (^-function has maximal type because ({s) has zeros at the negative even 
integers and 2^Tlog(T/27re) + O(logT) zeros with imaginary part < t < T and real part 
< (T < 1, so it fails both tests in Lindelof 's theorem. 

The Riemann ^-function is given by 

(2.2.4) ^(,) = i,(,_l)^-|T(i.)C(.). 

It satisfies the functional equation ^(s) = ^(1 — s). The ^-function is entire of order 1 and has 
the same non-real zeros as the (^-function, so it has maximal type. The Riemann Hypothesis 
is the assertion that all zeros of the ^-function have real part ^. 

Sometimes it is convenient to use the Riemann S-function, given by S(2;) = ^(| + iz). The 
S-function is an even function which is real on the real axis, is entire of order 1 and maximal 
type. Assuming the Riemann Hypothesis, S has only real zeros and n+(r) = n_{r) = 
(r/27r) log(r/27re) + O(logr). 

In the generalization of Theorem 12. 1.11 to entire functions of order 1 we will make use of 
the Hadamard product. Suppose / is entire of order p and let = [p] be the largest integer 
less then or equal to p. Then there exists a polynomial Q of degree at most A^ such that 

(2.2.5) f{z) = e«(^)z"° n (l - f ) (^^ (f ) ) ' 
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where Pn{z) = z + ^z'^ + ■ ■ ■ + j^z^ and Hq is the order of the zero of f{z) at z = 0. 

2.3. Generalization of Theorem 12.1.11 to entire functions of order 1. We have 

Theorem 2.3.1. Suppose f is an entire function of order 1 which is real on the real axis 
and has only real zeros. If q < p are consecutive zeros of f + af , a G M, then 

(2.3.1) inf Zn+i - Zn<p-q < sup Zn+i - Zn- 

If the zeros of f are simple and equality holds for one of the inequalities, then it also holds 
for the other, and f{z) = Ae^^ cos{Cz + D) for some real A, B, C , and D. 

So, if the zeros of / are not equally spaced, then differentiation makes the smallest gaps 
larger and the largest gaps smaller. It might not be obvious that if the zeros of / are equally 
spaced then the same holds for /', but one can check that the set of functions of the form 
Ae^^ cos{Cz + D) is closed under differentiation. In particular, the factor e^^ shifts all of 
the zeros of the derivative of cos(Ca; + D) by the same amount. 

Proof. We may assume that /(O) 7^ and the zeros of / are simple. We have the Hadamard 
factorization 

(2.3.2) f{z) = Ae^^ H " ^) {£) ' 
so 

(2.3.3) iz) = a + B + J2— • 

/ n ^ 

Imitating the previous proof, the terms a, B, and Zn cancel and equation ()2.1.4|1 becomes 

(2.3.4) = y i —. 

By (j2.1.3|) every term in the sum is nonnegative, and we have a contradiction unless p — q = 
Zj-^-i — Zj for all j. This proves the left hand inequality in the Theorem. The proof of the 
other inequality is identical. 

Finally, if p — g = Zj+i — Zj for all j then the zeros of / are equally spaced, so we recognize 
the Hadamard product of / to be of the stated form. □ 

A shortcoming of Theorem 12.3.11 is that the inequalities can be vacuous. For example, for 
the Riemann S-function inf Zn+i — Zn = 0- To formulate a "local" version of Theorem 12.3.11 
we must be able to speak of the "local density" of zeros. Suppose 

(2.3.5) n±(r) = L(r) + E(r), 

where L{r) is a nice function and E{r) = o{L{r)) should be thought of as an error term. 
Then L'{r) is the density of zeros near r, meaning that L'{r)~^ is the average gap between 
zeros near r. We want to show that if /' has zeros which are closer together than the average 
spacing, then in some neighborhood / will also have zeros which are closer together than 
average. The following Theorem asserts that the neighborhood can be taken to be of size 
E{r)/L'{r). 
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Theorem 2.3.2. Suppose f is an entire function of order 1 which is real on the real axis, 
has only real zeros, and n±{r) satisfies \2. 'J. 5]) where L is an increasing smooth function such 
that L^^\r) <^ r^^L(r). Suppose k < K < 1 and q < p are consecutive zeros of f + af , 
a e R, such that p — q < nL'^q)^^. Then there exist zeros z^, Zn+i of f such that \zn — q\ 
E{q)/L'{q) and Zn+i - Zn < KL'{qY^. 

For the Riemann S-function we have L{r) = (27r)^^r log(r/27re) and E{r) = logr, giving 
this 

Corollary 2.3.3. Suppose < k < K < 1 and 7^ < 72 are consecutive zeros of S*^") with 
72 — 7i < 27r/t/log7i. Then there exist consecutive zeros 71 < 72 of E with |7i — 7i| -C 1 
such that 72 — 7i < 2ttK/ log 7. 

The Corollary is motivated by the problem of Landau-Siegel zeros. Conrey and Iwaniec |3] 
have shown that if there exists < | and 6 > such that 72—71 < 2ttk,/ log 7 for ^ T log*^ T 
pairs of zeros < 72, 71 < T of the S-function, then one can find an effective lower bound 
for the class number of the imaginary quadratic field Q{\/—d). By the Corollary, if we had 
3> Tlog^ T pairs of zeros 72 — 7i < 27rK/ log 7 for for any fixed n > 0, then we would 
have ^ T such pairs of zeros for S. This falls short of the requirements for Conrey and 
Iwaniec's result. However, it may be possible to modify the proof of Theorem 12.3.21 to show 
that either S has a pair of zeros which are closely separated and very near to q, or else S 
has a large number of zeros which are not too far from q. 

Proof of Theorem \2.3.^A For a contradiction, suppose Zj^i — Zj > p — q ioi \q — Zj\ < X . We 
can reuse everything in the previous proof up to ()2.3.4|) : 







1 1 




(2.3.6) 
say. 

By assumption, 5*1 > 0. We now estimate S^, the treatment of 82 being identical. Note 
that 

1 1 1 1 P — Q 



P-Zj+i q-Zj q-Zj+i q-Zj {p - Zj+i){q - Zj+i) 
(2.3.7) = ^'\, +0^ ^""'^^ 



q-Zj+i q-Zj (g-Zj+i)2 \{q - Zj+if ^ 

Let Zj denote the smallest zero of / in [g + X, 00), so the sum in 6*3 begins with zj. In the 
sum 5*3, the first two terms in ()2.3.7|) telescope, giving 



(2.3.8) S2 = ^ + {p-q){l + 0{L'{qr) Y: . I 

Writing the sum as a Steiltjes integral and integrating by parts we have 



F — 

(n — 7 



Zj—q>X 



(g - Zj+iY J^j^^ (r - g)2 
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L'(r) , r dEir] 
^ ■dr+ ' ^ ' 



1 (r - qf J,j^, (r - q) 



2 



- q {zj+i - qY Jzj^, r-q J^^^^ (r - g)3 

(2-3-9) = (l + o(l))- - + - 

^j+i - q [zj+i - q) 

In the last step we used E{r) = o(L(r)) and L"(r) ^ L'{r)/r. 

Adding all the terms and using the fact that L'(zj+i) ~ L\q) and E{z,j^i) ~ -E(g) we 
have 

(2 3 10) ^3 ~ ^ _ (g-p)-^'(g) _ 

^ zj -q zj+i - q Zj+i - q 

We assumed q—p < KL'{q)^^, so the first two terms sum to something larger than 6/{zj+i—q) 
for some 6 > 0. The numerator of the last term is bounded by K,E{q)L'{q)~^X^^ which is 
smaller than 6 provided X > 6~^K,E{q)L'{q)~^. For such X we have 5*3 > 0, which gives a 
contradiction. We have X <^ E{q)/L'{q), as claimed. □ 

2.4. Long-term behavior. As described in the previous section, differentiation causes the 
small gaps between zeros of / to become larger and the large gaps to become smaller. Thus, 
one might expect that under repeated differentiation the zeros of Z*^"-' will approach equal 
spacing, which would equal the average spacing between zeros of /. And if the zeros of /^"^ 
are approaching equal spacing then f^^\z) should be approaching Ae^^ cos{Cz + D), where 
A, B, C, and D may depend on n. 

In order for the above argument to hold, there must be an appropriate sense of average 
spacing between zeros of /. Below we show that it is sufficient to have n+(r) ~ n_(r) ~ nr for 
some K > 0, but it is possible that this condition can be weakened. Since Ae^^ cos{Cz + D) 
has finite type, one might expect that this would also be a necessary condition on /, but this 
is not the case. The assumptions in the following result are stronger than the first condition 
in Lindelof's theorem, but weaker than the second. 

Theorem 2.4.1. Suppose f is an entire function of order 1 which is real on the real axis, 
has only real zeros, and ?T,+ (r) ~ n_(r) ~ nr. Then there exist sequences {An), (Bn), and 
{Dn) with Dn bounded, such that 

(2.4.1) lim Ane^"'f''''\K~^z + D„) = cos(7r2), 

ra— >oo 

uniformly for \z\ < X for any fixed X > 0. In particular, the zeros of approach equal 
spacing. 

We give the proof in Section |31 

It seems that An in the theorem is increasing or decreasing according to whether n{r) — 2nr 
is generally negative or positive, but we have not succeeded in finding a precise statement. 
The Theorem also suggests that the regularity in the distribution of the zeros of / should 
lead to a great regularity in the Taylor series coefficients of /. This seems worth exploring 
further. 

A trivial example of the theorem is when f{z) = Ae^^ cos{z + D) for some A, B, D. We 
have f{z) = VWTTf{z + <^) where cos{<^) = B/y/WTT and sin((^) = 1/VWTT. So the 
conclusion of the theorem holds with An = A^^{B'^ + 1)^2" and Dn = —D — mp (mod27r). 
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Somewhat more generally, suppose f{z) = Ae^^ (cos(2; + D) + E) where \E\ < 1. This 
function has zeros consisting of two interlaced sequences which separately are equally spaced. 
The conclusion of the theorem holds with the same An and D„ as in the previous example. 
Similarly, if f{z) = F{z) cos(z) where F is a degree N polynomial, then the conclusion holds 
with A„ = n-^ and G {-f , 0, f , vr}. 

A less trivial example is the Bessel function Jn{z). We have n±{r) = vr^^r + 0(1), so the 
theorem applies. In this case the theorem can be verified directly because 

1 r 

Jn{T^z) = — cos{nzsmt — nt) dt, 

Jo 

SO for an even number of derivatives 

jj^^'^^TTz) = / sin^'' (t) cos{n z sin t-nt)dt 

Jo 

(2.4.2) ~ (_^)^-^cos(7r2;-^), 

because sin tdt is approaching a Dirac 5-function at t = 7r/2. So the conclusion of 

the theorem holds with A2k = \/'kpK and = -y- mod 27r, and similarly for k odd. 
Equation ()2.4.2|) can also be obtained from the Taylor series for J„. 

Theorem 12.4.11 asserts that the zeros of /^"^ approach equal spacing. The approach to 
equal spacing is actually quite fast. The general case is somewhat intractable, so we will 
assume that the zeros are close to equally spaced and do a first order approximation. We 
also suppose further that / is approximately an odd function, although this is not essential 
and we describe the necessary modifications during the proof. Thus, 

(2.4.3) /(,) = C(.-.„)n(l--f^)exp(-^). 

where \ej\ < e for some small e. 

Theorem 2.4.2. Suppose f is given in \2.4-^ and let zf^ be the zeros of the jth derivative 
/(^■). Thenzi^,-zi = l + 0{er'). 

We give the proof in Section 14.31 

If the zeros of / do not have a nonzero average spacing, then the results in this section do 
not apply. Suppose, for example, n^{r) and n_{r) grow faster than linearly. The "overall" 
average spacing is 0, and under repeated differentiation there is a competition between zeros 
moving toward the origin and zeros trying to become locally equally spaced. If n+(r) and 
n_(r) are sufficiently nice then there is the possibility that an analogue of Theorem 12 . 4 . 1 1 mav 
hold. This was established by Haseo Ki [S] in the special case of the Riemann S-function: 

Theorem 2.4.3. (Haseo Ki ^) There exist sequences An and Cn, with ^ slowly, such 
that 

(2.4.4) lim An {Cnz) = cos(z) 

n— »oo 

uniformly on compact subsets of C. 

The above result was conjectured in an earlier verision of this paper. The S-function 
has n±{r) = (27r)^^r log(r/27re) + O(logr) so the zero counting function is particularly 
well-behaved. However, Ki's proof makes use of the fact that the S-function has a nice 
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representation as a Fourier transform, rather than directly using properties of the zeros. 
The theorem suggests that log(^^")(|)) should grow very regularly and not too much faster 
than linearly a.s n —>■ oo. Rick Kreminski [121 E] has calculated the first 490 derivatives and 
his data appears to grow approximately like nlog?7.. 
More generally we conjecture: 

Conjecture 2.4.4. Suppose f is a real entire function of order 1 having only real zeros, 
with n+(r) and n^{r) sufficiently nice and n^{r) — n_{r) not too large. Then there exist 
sequences An, Bn, Cn, and Dn with Dn/Cn hounded, such that 
(2.4.5) lim An e^"7(2n) ^^^^ ^ _ ^^g^^^, 

uniformly on compact subsets of C 

3. Averaging, instead of differentiating 

Differentiation is a process which takes one sequence of points (the zeros of /) and replaces 
it with another sequence of points (zeros of /'), such that the two sequences interlace each 
other. For comparison, we also consider the much simpler process of making a new sequence 
from the midpoints of neighboring elements of the given sequence. 

3.1. Averaging on the line. Suppose (x„) is an increasing sequence. Form new sequences 

{xD where x° = x„ and x{^^ = ^ + x-'^_^_(^_iy^ ■ That is, the terms of each new sequence 

are the midpoints of consecutive terms in the previous sequence. The subscript n + (—1)-' 
is designed so that if (x^) is equally spaced, then x^"*"^ = xi^. If that subscript was n + 1 
then the equal spaced case would give x{^^ = and the terms would be "drifting to the 
right." 

Theorem 3.1.1. Suppose {xn)n& is a sequence with Xn = n + En, where En <^ E{\n\), as 
n —>■ ±oo, for some increasing function E. If E{n) = o{n), that is Xn ~ n, then as j ^ oo 
the j^^ midpoint sequence (x{) approaches equal spacing. Furthermore, if E{n) <^ with 
0<9 <1 then 

<+i-< = l + 0((n^+j^')r^ 
Proof. It is straightforward to show by induction that 



J2 [l)^n+[j/2]-k- 



0<k<j 



Thus, 



^ ^ 5Z {^n+l+[j/2]-k - Xn+[j/2]-k) 

0<k<j ^ ^ 

Yl [i)^'^^ 5Z [Ij {^n+l+[j/2]-k- £n+lj/2]-k) 



0<A:<j ^ ^ 0<k<j 

'J\ (J 



2 ^ Y ^n+l+[j/2]-k ( K 



0<fe<j 



k \k~l 



+ 2 ^ {en+[j/2]~j - ^n+[j/2]) 



DIFFERENTIATION EVENS OUT ZERO SPACINGS 



11 



(3.1.1) = 1 + 2- ^"+4zlJr^7^ +^(2-(n + j)). 



2 ' 2 



We may assume that E is an increasing function and E{a + b) ^ E{a) + E{b) for a,b > 0. 
So the sum in ()3.1.H) is bounded by 

(3.1.2) J2 iE{n)+E{i)) (^^1^) 



|-£+l 



Now use fact that the binomial distribution approaches the Gaussian to estimate the 
binomial coefficient: 



(3.1.3) 

So, 




2^ / 2f 

exp 



0<^<ii ^ ' 1 J 0<£<i- ^ 



- 2 



< — / E(t) exp ; dt 



r Jo V J y |-t + l 
(3.1.4) = [ ' Eity^)exp{-2t^)- ^-^-dt. 



We may assume E{ab) <^ E{a)E{b) for a, 6 > 1, so ()3.1.4p is bounded by 

(3.1.5) 2^rE{^j) / exp(-2t2)^ dt 2^ E{^), 

Jo ^ - tv.^ + 1 

the integral being seen to be -C by breaking it at j^, for any < 5 < |, and estimate 
the two parts separately. 

To finish the proof, note that if E{n) = o{n) then j~^E{y^) = o(l), and if E{n) ^ 
then the other estimates follows immediately. □ 

If x„ = n + 0(1) then the above shows that the jth midpoint sequence approaches equal 
spacing with a discrepancy of order j~2. This convergence to equal spacing is actually 
quite slow, and in some sense it is the slowest way to "even out" a sequence. In particular. 
Theorem 12.4.21 savs taking successive derivatives evens out the sequence much faster, with a 
discrepancy of order These topics are discussed in Section lOl 

3.2. Averaging on the circle. As mentioned in the Introduction, it is believed that there 
is a close analogy between L-functions, which (conjecturally) have all their zeros on a line, 
and the characteristic polynomials of random unitary matrices, which have all their zeros 
on a circle. See '3]. It is obvious that the analogy must break down if pushed too far, for a 
polynomial has only finitely many zeros. In this section we consider the averaging of points 
on a circle, in analogy to the averaging on a line in the previous section. Here we find that 
the discrepancy from equal spacing vanishes exponentially, in sharp contrast to the j~2 of 
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the linear case. This exponential convergence is a general property of iterative averaging 
procedures on finite sets, as described in Section 4.4 of [7j. 

Suppose pi, . . . ,pi are points on the unit circle, and form new sets of points on the circle 
in analogy to the previous section, where < p:^ < 27r is interpreted as an angle. If j < n 
then this process is indistinguishable from the averaging process in the previous section. But 
when j > n the fact that there are only n points on a circle comes into play. This causes the 
averaging to even out the spacing more rapidly. 

Proposition 3.2.1. Let and be sequences of i points on the unit circle as described 
above. Then p>^^, - p^^ = ^ + 0((1 - ^Y). 

Proof. Starting exactly as before, but using the fact that pj = pk if j = k mod i, we have 

PL+I-P'n = Yl [I) {Pn+l+[j/2]~k-Pn+[j/2]-k) 



0<fc<i 

e 



2 ^ ^iPm+l -Prn) 

m=l 0<k<j 

k = M mod e 



where M = n + 1 + [j /2] — m. 

Now let p be a primitive i^^ root of 1, and let pj = />' for < j < i — Using the binomial 
theorem and the fact that Pfc = except when n = mod i, we have 

(3.2.1) i Yl (^) = ii + iy+PM,a+piy + ---+PM,.,a+Pi-iy 



0<k<j 
k = M mod £ 



(3.2.2) = 2^+0 2^£1- — 




where Pm^ is an ordering of the pk depending on M. The last step follows from |1 + e^^^^^] 



2(1-^) +0(1- 

Finally, use the fact that J2m=i (Pm+i — Pm) = 27r to finish the proof. □ 

4. Almost equally spaced zeros 

Suppose f{z) is an entire function of order 1 which is real on the real axis, has only real 
zeros, and the zeros have average spacing 1. As we repeatedly differentiate, the zeros will 
approach equal spacing, and we want to determine the rate at which this occurs. We will see 
that the approach to equal spacing is much faster than in the midpoint process of Section El 

We treat in detail the case that the zeros of / are close to equally spaced, both on the line 
and on the circle. At the end of this section we determine the rate at which differentiation 
evens out zero spacing. 

4.1. Almost equally spaced zeros on a line. Suppose 

(4.1.1) „.).C,.-.„,n(l--^)exp(-^). 

where \ej\ < e for some small e. That is, / is approximately an odd function with zeros close 
to equally spaced. 



DIFFERENTIATION EVENS OUT ZERO SPACINGS 13 

Theorem 4.1.1. In the notation above, f has zeros at k + ^ + a^, where 
and f" has zeros at k + (3k, where 

3 TT^^ j2 

In the proof of Theorem 14.1.11 we wiU require the foUowing standard formulas: 



(4.1.2) y"- 777 = —, and V 

^ ' ^ ill + 124' ^ 



8?i2 

^ (2j + l)2 4' ""^ Z^(2j + l)2(2(j+n) + l)2 ^ = 0- 



Proof of Theorem We have 



f 



Here and below, means that the j = term must be modified in the obvious way. The 
zeros of f'{z) will be approximately halfway between the zeros of /, so suppose one of the 
zeros is k + ^ + ak- 

Note: if we didn't assume that / was odd, and so a factor of e"^ occurs in our function, 
then the zeros of /' would not be halfway between the zeros of /. But the zeros of /' would 
all be shifted from the midpoint by the same amount, and all of the calculations below would 
work with a slight modification. So, 

= y(fc + ^ + «fc) = E',^ , 

1 + 
4 

£j+k 



f 2 ^ + I + ftfe - J - e:^ j + ej 

^l + 2(k- 1) ^ (1 + 2(k - iW 7 + ^ ^ ^ ' 
J2 n^'^'- ^^.^^. +0(6' + at 



(1 + 2(A;-j))2 



(4.1.3) = 4^5^ ^27^-"'"^ + ^^"' + "^)- 

Solving for ak gives the first formula in the Theorem. The 0(^2) error term follows because 
the above formula shows that ak = 0{e + e"^ + al), but = o{e) and al = o{ak), so 
ak = 0{e). 

To understand the effect of the second derivative on the zeros, we must iterate the above 
formula. The second derivative /" will have zeros near the integers, so suppose there is a 
zero at + jSk- By the first formula in the Theorem, and using formulas ()4.1.2|1 . 

4 aj+k^i 



7r2 ^ (2j - 1)2 7r2 ^ (2m - 1)2 ^ ^ ^ 
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16 ^ 1 1 2 

—A 1^ ^n+k 2^ (2j_l)2(2(n-J + l)-l)2 + 



3 

(4.1.4) = 1 +^^£^ + 0(a 

as claimed. 



□ 



4.2. Almost equally spaced zeros on a circle. If the polynomial / has all its zeros on 
a circle, then /' has all its zeros strictly inside the circle, except at points where / has a 
multiple zero. This follows from the famous Gauss-Lucas theorem that the zeros of /' lie 
inside the convex hull of the zeros of /. Thus, if we are to find an analogy with the case of 
zeros on a line, we must do something slightly different than differentiation. 

Lemma 4.2.1. If f is a degree n polynomial with all its zeros on the unit circle, then 

also has all its zeros on the unit circle. The zeros of f^~^^ interlace the zeros of f^ . 
Proof. From 

(4.2.1) f{z)= J2 «fc-2' = «n n ('^-^^■) 

0<fc<n 0<A:<?i 



and the fact that l/zj = 'zj, we have f{z) = (ao/a„)2;"/ (1/^). It follows that Z{z) : = 
{ao/an)~^ f(z) is real on the unit circle. By RoUe's theorem, ^Z(e*^) has n zeros in 
< 6 < 27T. But these give the n zeros of Z', which are the same as the zeros of □ 

We now derive the analogue of Theorem 14.1.11 Suppose the zeros of / are close to equally 
spaced on the unit circle: 



where \ej\ < e and we set e(x) = e^'^*^. In the above formula and in Proposition I4.2.2| 
unrestricted sums and products over j should be interpreted as over j modulo n, and we 
consider Ej to only depend on j modulo n. 
Let 

g[z) = Z-'^f{z) 

and suppose g' has a zero aX e{^ + ^ + au). 
Proposition 4.2.2. In the notation above, 

Note that the above sum is 

, oo 

Sj+k 



E 



7r2 ^ {23 -ly 
as n — >■ oo, so this formula matches that in Theorem 14.1.11 
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^2 



Corollary 4.2.3. C(2) 
Proof. Using sin(x)~^ = + 0(1), from Proposition I4.2.21 we have 



7r2 ^ (2j - 1)2 ' I n2 

4 V-^ 



2 ^J- 2 



as n — i> oo and e 0. 

If Ej = e for all j , then by a change of variables we see that the zeros of g' are ate(^ + ^+5). 
This means that at = e, so canceling e in ()4.2.2|) gives 



= 2 

(2j~i)2 ' "which is equivaicui uu z^j=i ^ 



Since n and e are arbitrary, 1 = (2^1)^' "which is equivalent to Yl'jLi ji = D 



Proof of Proposition \4-2.^ We have 



^(-) = -- + y 



so setting 2; = + ^ + a^) gives 



n 2n 

n 



Write -E'(t) = e{t) — 1, so that it is easier to keep track of the main terms: 



, ^(^ \E{e,)\' + \E{a,)\' 
The first terms on both sides of the above can be canceled, because both equal 



dz ^ ' 2" - 1 ^ 



evaluated at e(- + 77-). 

Rearranging to put the terms with on one side and multiplying by e(^ + 2^), we have 



n 
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The sum on the left side equals _!i(lL-2) because it equals 



-z'—iog{z- - 1) = — = E 



dz^ ' (z^-iy ^(^-e(l))2 

evaluated at e(- + 77-). This brings us to 

-E(-ak) - -(^2 - 2n)E(ak) = Y r , . ^^^'^ , , . + O (nV) 

j v'^V2n ^ An 2n> 2n An ^ 2n)> 

Now use the approximation E{t) = t + Oif'), for t small, and solve for a^. 



Oik 



1 V + 0(r 

,2 Z^™2/^2i-lx 



sm , „ 

3 ^ 2 71 



as claimed. □ 

4.3. Rate of convergence to equal spacing. Now we analyze the rate at which differen- 
tiation evens out zero spacings. 

For both the midpoint and differentiation process, applying the operation twice has some- 
what nicer properties than just applying it once. For example, if the points are equally 
spaced then applying the operation twice leaves the sequence unchanged. In the discus- 
sion below we will generally work with the second derivative and second midpoint process, 
although the same ideas apply to the basic process. 

We will view differentiation and midpoint averaging as two examples of a general proce- 
dure. Suppose P is a probability measure on the integers such that P is even, P(0) is the 
maximum, and P{j) is increasing for j < and decreasing for j > 0. The conditions on 
P are natural if one is thinking in terms of using P to smooth out the irregularities of a 
sequence. 

Given a sequence (e^) define a sequence of sequences by /5° = e„ and 



(4.3.1) = E^H/^" 



For the second derivative process -P(O) = | and P{n) = -^^^^^ otherwise, and for the (second) 
midpoint averaging process -P(O) = |, P{—1) = -P(l) = \, and P{n) = otherwise. 

In Theorem 12.4.21 we start with Zk = k + Ek, so in the notation above the discrepancy of 
z^_^_-^ and z^ from the average spacing is — Pi- We wish to estimate this in terms of £„. 
We can write Pi in terms of e„ as 

(4.3.2) (3i = J2P\n)en+k, 
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say. We can estimate the difference as 

f^Ui-f^i = Y.'n+k{P\n)-P\n-l)) 

n 

< £^|P^(n)-P^(n-l)| 

n 

(4.3.3) = 2eP^{0). 

The last step requires that P^(0) is the maximum, and P^{n) is increasing for n < and 
decreasing for n > 0. We will see that this is the case. 

It can be seen that P^ = P*P*---*P is just the £-fold iterated convolution of P with 
itself, where 

(4.3.4) P*Q{n) = J2Pi'^)Qi^-^)- 

m 

The 1 = 2 case was demonstrated in Section f4.11 where was calculated from a^. It is 
straightforward to check that if is the global maximum for P and Q, and both functions 
are increasing for < and decreasing for > 0, then the same holds for P * Q. Thus, 
inequahty ()4.3.3|1 applies and we need only evaluate P^(0) = P * P * ■ ■ ■ * P(0). 

The evaluation can be done using the properties of convolutions and Fourier transforms. 
Given a sequence S, let 

(4.3.5) Fs{x) = 5(^)e^™^ 

n 

Then 

Fs{x)FTix) = (^X^5(n)e2™^j (^nj)e'-^^^ 

{n+j)x 



n j 

= J2e^~J2S{n)T{m- 

m n 

(4.3.6) = ^(5*T)(m)e2~. 



That is, convolution of sequences corresponds to multiplication of Fourier series. 
The final ingredient is to note that if 

P(a:) = J]c,e2-^^ 
i 

then ^ 

Co = / F{x) dx, 
Jo 

provided we can integrate term-by-term, because e^^^^^ dx = unless j = 0. Thus, 

P^(0) = P*P*---*P(0) 

(4.3.7) = [ (^P(n)e2™^) dx. 
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Now we need only identify the function in the integrand above. 
In the case of the (second) midpoint process, we have 

J2 P(ri)e2™^ ^l^l cos(27rx). 

n 

You can either recognize the Beta integral or ask a computer algebra package to verify that 

1 ^1 1 



^ (0) = / (2 + 2 ^^^^^''^^ ) 

2r(f + £) 



(4.3.8) 



;i + 2£) v/?r(i 
1 



0(f 



which is the error term from Theorem 13.1.11 assuming -C 1. 

In the case of the (second) differentiation process, it follows from the Fourier expansion 

/ 1 4 ^ cos(27rnx) ^ ^ 

^ ' n=l 



that 



for < X < 1. Thus 



2 



^P(n)e^— = 4('x-^ 

n ^ 



P'{0) = / 4 rfa; 




(4.3.9) 



as claimed in Theorem 12.4.21 

The discrepancy from equal spacing is much smaller for the differentiation process than for 
the midpoint process, and this is not surprising. Both processes average among neighboring 
points, and that averaging will be more effective if takes place over a larger range. In the 
setup we have described, if P has finite variance cx^ then Fp{x) ~ 1 — 27r^(T^x^ for x near 0. 
Therefore 

' 1 



Fp{x) dx 



^5 V2TTaH 

as we saw for the midpoint process. For the differentiation process we found that Fp{x) 
1 — b\x\ + 0(a;^) for x near 0, where 6 = 4. Therefore 



^ 2 
Fp{xYdx^ — 

-s 



bt 
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as in Theorem 12.4.21 Many other behaviors are possible. If P has infinite variance but 
decreases hke then Fp[x) ~ 1 — a|a;|" as x — * 0, where a = A — 1, and so 




Thus, we have averaging processes where the £th iterate approaches equal spacing as fast as 
any given power of i. 

It should be noted that if one violates the conditions that -P(O) is a maximum and P{n) 
increases (decreases) for n < (n > 0), then the "averaging process" may not lead to 
equal spacing. For example, if -P(— 1) = -P(l) = | then the process leads to two interlaced 
sequences which separately approach equal spacing. 

5. Proof of Theorem 12.4.11 

We wish to show that if / is an entire function that meets suitable additional conditions, 
then the nth derivative /^"\ appropriately rescaled, approaches the cosine function. The 
real issue here is proving that the zeros approach equal spacing with sufficient uniformity. 
For if the zeros approach equal spacing, then the Hadamard product for /^"^ can be seen to 
be close to Ae^^ cos(Cx + D) for some A, B, C, D. So we first show that the zeros of f^'^^ 
approach equal spacing, and then the proof is almost immediate. 

Our method makes extensive use of the midpoint process we studied in this paper. It 
would be interesting to find a more direct proof, which could possibly lead to a stronger 
result. At the end of this section we also discuss the case that the zeros of / lie near, but 
not necessarily on, a line. 

5.1. Repeated differentiation leads to equal spaced zeros. We have shown that under 
repeated differentiation the small gaps between zeros are becoming larger, and the large 
gaps are becoming smaller, but it does not trivially follow that those gaps are approximately 
equal. Since we know that the midpoint process gives gaps which approach equal spacing, 
one possibility is to show that differentiation is better than midpoint at evening out the 
sequence: 

Conjecture 5.1.1. Differentiation is better than midpoint. Suppose f is an entire 
function of order 1 which is real on the real axis and has only real zeros, and suppose Zj are 
the zeros of f , listed in increasing order. If q < p are consecutive zeros of f then 

1 1 

inf -{Zn+2 - Zn) <p-q < sup -{Zn+2 - Zn). 

A similar result should also hold when the spacing between zeros of / varies slowly, in 
analogy to the relation between Theorem 12.3.11 and Theorem 12.3.21 Note that we interpret 
the right side of the above inequality as "oo" if / has only finitely many zeros. 

In contrast to Theorems 12 . 3 . 1 1 and 12 . 3 . 2l the conjecture is not true if instead we assume p 
and q are zeros of /' + a/, for if a is large then p and q are very close to zeros of /. Also, 
the conjecture is not true if |(2n+2 — Zn) is replaced by |(z„+3 — Zn)- 

We now show that repeated differentiation leads to equal spaced zeros. The proof makes 
use of the fact that the midpoint process leads to equal spacing, but our approach is somewhat 
less elegant than would follow from the above conjecture. 
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Proposition 5.1.2. Suppose f is an entire function which is real on the real line, has only 
real zeros, and n+(r) ~ n_{r) ~ r. Then if (z^) are the ordered zeros of the jth derivative 
f^^'^ then zl^j^i — z^ = 1 + o(l) as j oo. In addition, \z-^\ ^ n. 

Proof. The same method as in the proof of Theorems 12.3.11 and 12.3.21 shows that if {zn) are 
the ordered zeros of / and (z^) are the ordered zeros of /', then 

(5.1.1) inf(2;„+Ar - Zn) < - 4 < SUp(z„+Ar - Zn), 

for any A^, and similarly when considering zeros with \zn — z'^ < X. When N = 2, the 
difference Zn+2 — is twice the gap between the midpoints of consecutive zeros of /. By 
Theorem 13. l.H iterating the midpoint process gives sequences which approach equal spacing. 
Thus, the sequence of next nearest neighbors of zeros of f^^^ approaches equal spacing. That 
is, the zeros of f^^^ consist of two interlaced sequences which separately are approaching the 
(same) equal spacing. 

It remains to show that the whole sequence of f^^^ zeros is approaching equal spacing. One 
possibility is to now prove that /(^^ {z) is approximately of the form Ae^^ {cos{Cz + D) + E), 
so then the zeros of derivatives of f^^^ will approach equal spacing, as discussed in the 
paragraph after Theorem 12.4.11 But our goal is to prove that the zeros become equally 
spaced and then use that as a tool to prove the functional form. 

We will show that the zeros of f^^^ also consist of three interlaced sequences, each separately 
approaching equal spacing. From this it immediately follows that the entire sequence is 
approaching equal spacing. We will use (jS.l.ip with = 3. Note that 

Zn+3 ~ Zn = {Zn+3 + + ^n+l) "~ (^n+2 + + ^n) 

(5.1.2) = 3{Zn+2 — Zn+l), 

where Zn = |(-Zn-i + Zn + Zn+i)- By either imitating the proof of Theorem 13. 1.11 or using the 
method of Section lOl with -P(— 1) = -P(O) = -P(l) = |, we see that iterating the averaging 
process (zn) (5„) leads the sequence to approach equal spacing. So zl^_^_^—zl is approaching 
equal spacing, as required. 

The final assertion follows from the fact that z° ~ n and that each successive zero set 
interlaces the previous one. □ 

5.2. Proof of Theorem 12.4. 11 Assume k = 1, so the zeros of f^^^ are 1 apart on average. 
We first shift the function slightly so that the zeros are arranged conveniently. Suppose the 
smallest non-negative zero of /'••'^ is at wi and the largest negative zero is at Choose 
dj = |(w_i+Wi) and let zi, Z2,- ■ ■ denote the positive zeros of f^^\z + dj) and z^i, z„2,. . .the 
negative zeros. 

We have the Hadamard factorization 

f^^\z + d,) = A,exp{B,z)l[(l~-)exp(- 

V Zn I \ Zfi 

n ^ / \ 

= Aexp(5;.) n (i-;^)(i-f 



\<n<Y 
X 



n (i-f )-p 

\n\>Y 



Z 

Zn. I \ Zn. 



(5.2.1) = A,^x^{B]z)F,{z)G^\ 



z 
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say, where Y is to be chosen later. We must show that if \z\ < X then Fj{z) cos(7r2) and 
Gj{z) 1. 

Given X, e > 0, first choose Y so that \Gj{z) — 1| < ee~'^^ for \z\ < X and all j > 0. 
This is possible because Zn ^ n. 

Let F and G denote the above products in the special case that z±n = -F |)~^ for all 
n > 0. Since F{z)G{z) = cos{nz), it follows that | cos(7rz) — F{z)\ < e for \z\ < X. 

For each n, as j — » oo the gap 2;^ — is approaching 1, so we can a choose j so that 
z±ri is sufficiently close to {n =f |)~^ for < n < F to ensure that \Fj{z) — F{z)\ < e 
for \z\ < X. Using |cos(7r2;)| < e^'^ for \z\ < X and combining all estimates we have 
\A~^ ex\){—B'-z) f'^^^z + Dj) — cos(7r2;)| < Ae for \z\ < X, which completes the proof. 

5.3. Zeros not on a line. We have seen that repeated differentiation leads to a function 
whose zeros approach equal spacing. Surprisingly, this does not always require the zeros 
of the original function to lie on a line (or circle), but only that the zeros lie in a suitable 
neighborhood of the line (or circle). 

This is easy to see in the case of the circle. Suppose / is a degree n polynomial and let 
g{z) = z-^Z-'fiz), so 

(5.3.1) g{z) = anZ2 + an-iZ2 + ■ ■ ■ uqZ 2. 
Then 

(5.3.2) zg {z) = -anZ2 + an-iZ2 -\ \^--jaoz ^, 

and the "fcth derivative" equals 

(5.3.3) (^^) ^^^^^[2) «"^' + (^2~V V 2/ 

If Oq 7^ then the factors (1)*^ from the first and last terms dominate everything else, so the 
zeros are approaching the zeros of + (—1)*^^. Note that we didn't actually require the 
zeros of the original polynomial / to be on the unit circle. As long as is not a root of /, 
the above process gives a sequence of functions whose zeros eventually lie on a circle, and in 
fact approach equal spacing on a circle! 

For zeros on a line, we have the following result of Young-One Kim (this is a slight revision 
of Theorem 2 from ^2] ) '■ 

Theorem 5.3.1. Let f{z) be a nonconstant real entire function, < p < 2, and assume 
that f{z) is of order less than p or is of order p and minimal type. If there is a positive real 
number A such that all the zeros of f{z) are distributed in the infinite strip \Qz\ < A, then 
for any positive constant B there is a positive integer ni such that f^'^\z) has only real zeros 
in |3fJz| < Bn~p for all n >ni. 

In other words, if the zeros don't start out too far from the real axis, then they end up 
on the real axis as you differentiate. Thus, Theorem 12.4.11 can be modified to only assume 
that the zeros lie in a strip around the real axis, and Conjecture 12.4.41 should only require 
the zeros to lie near the real axis. 

The theorem applies to the Riemann S-function with p = 1 + e. It is interesting to note 
that Conrey has shown that S*^-^^ has a positive proportion of its zeros on the real axis, 
and that proportion is 1 + 0{j^^) as j — > 00. 
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If f{z) = Xlj^o '^i-^Vi' is entire of order 1 and all Cj are real, then a necessary condition 
for the zeros of / to be real is the Turan inequalities — Ck-iCk+i > 0, for A; > 1. The above 
ideas suggest that the Turan inequalities should hold for sufficiently large k, provided only 
that the zeros of / lie in a neighborhood of the real axis. 
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